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On Interpolation, Summation, and the Adjustment of Numerical 
Tables. By W. S. B. Woolhoxtse, F.R.A.S. ; F.I.A., &c. 

Summation. 
1 HE three heads, into which it has been found expedient to divide 
the general subject of the present paper, all relate to the consideration 
of a series of numerical quantities which are supposed to represent, 
either accurately or approximately, some known or implied law 
which necessarily has an influence in regulating their progression. 
Assuming that we only know the values of certain terms of the 
series of quantities, distributed at stated equidistant intervals, it 
has been shown that all the intermediate or absent terms are most 
accurately and expeditiously determined by the process of inter- 
polation, already discussed and practically developed. In some 
calculations, however, the intermediate values may be considered 
to be mere auxiliaries, and the ulterior object of the interpolation 
may simply be that of obtaining the sum of the complete series of 
quantities so found. Now, the final result of this last summation 
may be directly determined from the given values and their differ- 
ences by a special formula, without incurring the labour of inter- 
polating the absent terms of the series. It is this particular 
formula, and some remarkable and useful modifications of it, which 
we now propose to investigate. 

The general expression (j) for an interpolated quantity V, is 
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V, = {(V) + (a:-i)„ 1 } + ±I± ) {(S) + l=i Cl | 

+ 27^6 l (/)+ ~l~ 9i J 

, T (, 2 -l X ^-4)(^-9X,-4) i *__! , &c _ 

2 8 l v 7 9 J 

On examination, it will here be found that the several coefficients 
have like values when x becomes 1 —x; but those of a x , c lt e l} g x , 
»!, will appear with reverse signs. Thus, for example, the coefficient 
of (/i may be written 

( x -3)( x —2)(x-l)*(x+l)(x + 2) x— \ 
2. ...6 ■ 7 ' 

and the several factors obviously recur with a contrary algebraic 
sign when x is replaced by 1 — x. It is, therefore, evident that the 
result of either a summation or an integral, from x=0 to#=l, 
will vanish with respect to «,, c„ e { , g x , «',, the coefficients of which 
contain an odd number of these changing factors. 

Let m be the number of subdivisions of the interpolation, so 
that m— 1 values of V are interposed between V and V,, and the 
number of values will be m+1 if V and V, be included. The 
corresponding m+1 values of x are 

12 3 to 

u, — i — , , . . . . . 

to to TO TO 

To proceed with the summation, it will first be requisite to 
determine expressions for the sums of the powers of the numbers 
1, 2, 3, .... m. Let the prefix 2 denote a summation of the 
values of any quantity for all integral values of m from 1 to m, or 
from m down to 1. Then, in accordance with this definition, 
2{, w »+i _( m _ i)»+i } =OT "+i —(to- 1 y+i 

+ (■»»— 1)"+i_(to— 2) n+1 

+ (to— 2)"+'— &c. 
&c. &c. 

+ 2"+>-l 

+ 1-0 

In place of the quantity under 2, substitute the equivalent 
expansion 

(w+i V- ' to--'+ v ^ / ; - ; to- 2 - &c, 

and divide by n+ 1, then 



1864] On Interpolation, Summation, fyc. 303 

*»■- J JW- + =&=!> Sm- 2 -&c.= =£ 
^ *«o w-p J. 

.-. Sw"= ; + -Sot— 1 ' - ' Sw— 2 

»+l 2 2.3 

»(«— 1)(»— 2)„ „ . . 
+ V 23 A 4 y 2m*-«-&c. 

where the factors after the first term are the binomial coefficients 
of the order n, divided respectively by 2, 3, 4, &c. 

Hence, putting 1, 2, 3, &c., successively for n, we obtain 



*Sm= 



Zm 2 = 



2 
2w 3 + 3m 2 +»n 



Sm 4 = 

Sm 5 = 



6 
m' l -\-2m z +m? 



4 
6m 5 + 15m 4 + 10m 3 — m 

30 
2m 6 +6m 5 +5»* 4 -m 2 



12 

fl _ 6w 7 +21>n 6 + 2lOT 5 — 7m 3 +m 

„ , 3m 8 + 12m 7 +14m 6 — 7m 4 +2m? 
Sw'= 



2m 8 = 



24 

10m 9 +45m 8 +60m 7 — 427w 5 +20m3-3m 
90 



Ifi= 


=2m + l; 


4=m(m + l); then 










Sm = 


2 










£m 2 = 


kk 
6 










Sm» = 


fc 2 










Sm 4 = 


'«<»- 


i) 








Sm 6 = 




i) 








Sm 6 = 


=S<»- 


-34 


+ 1) 



2m 7 = |^(34 2 -44 + 2) 
2m 8 = ^ (54 s - 10J 2 + 94 - 3) 

2m 9 = |jr (2/t 3 - SF + 64 - 3). 

The following properties are readily inferred : — 
(2m) 2 =SW 
(Sm 3 ) 2 = J(2i» 5 + SmO. 
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j" V 



Now, for the sake of symmetry and the facility of carrying on 
the summation into the adjoining intervals, it will be desirable to 
reject one-half of the extreme values V , V,, for which *=0, 1, 
respectively. 

For example, suppose that m=5, 
and in the diagram let the ordinates 
PM, QN, represent geometrically the 
values of V , V,, and pm, p'm',p"m", 
p'"m'" the interpolated values denoted 
by v, v', v", v'"., Then the summa- 
tion under the proposed modification 
will be 

2(VV.. 1 =iV +f+t/+»"+t..'" + !V l 

I I II II , 

r -\-v v + v 




M 



V 



jr* 



v, 



W W m* - ^ 



2 + 2 



+ 



r'" + V, 



each term of which now expresses a mean value (V) of two con- 
secutive ordinates. 

Let the extremities of the ordinates be joined consecutively, so 
as to form the polygonal line P p p' p" p'" Q, and let the trapezoidal 
areas contained between the successive ordinates be denoted by 
oh °2j «3> «4> «»• Then, assuming that Mm=mm'=&c. = l, we 
shall have 



V + r 



=a,, 



v + v 



■ =«2> 



=a 3 , &c; 



and hence 



2 C v ")o....l=ai + "2+< l 3 + a4+a5=P'- 



Thus, supposing the distance between the interpolated ordinates 
to be the linear unit, we conclude generally that when, as proposed, 
one-half of the extreme terms are rejected from the summation, the 
result will express, geometrically, the polygonal area contained 
between V and V, . 

By modifying accordingly the signification of S(x n ), that is, by 
including only one-half of the extreme terms in the summation, 
and distinguishing the same by means of a parenthesis, observing 
that the values of x involved in the summation are 



o, 


1 


2 


3 


m 


m 


™' 


» * • 
m 


m 



we shall then obviously have 



s ( , )= £- - j( „. +1 , = --_ 1; 
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and thus, from the preceding expressions (4), we deduce 
2(l)=m 
. . m 

S(*«)= 



5 T 3m 30m 8 

. .. m 5 1 

*■ - >— 6 + 12m "" 12m 8 

. „_ ml 1 1 
W— 7 + 2w" - 6ot» + 42^ 

„/ « »» 7 7 1 
2(* 7 )=7T + 



8 12m 24m 8 ' 12m s 
m 2 7 2 



S (^)=o+5IT-T7Z3 + 



9 T 3m 15m 8T 9m s 30m 7 ' 

We now proceed to apply these to the foregoing formula for 
V x , in which it has already been explained that the result will be 
evanescent as regards a i} c 1} e 1; &c. 

The other terms, when expanded, are 

(V)+^-V-*) 

< (^ ;(««— S* 8 — 5«*+15« 8 +4a;»— 12*) 



2... .6' 



+ -i^I— (a8-4a: 7 - 14^+56^ + 49* 4 — 196x8-36^+ 144x). 
2 . . * .8 

The severally corresponding coefficients for the required sum- 
mation are therefore deduced as follows : — 

Coefficient of (V)=m. 



•itto of ^ = 


m 

3 

m 

~2 


+ 


1 
6m 




m 




1 


— 


~ 6 


+ 


6m 



n'—l 
6m 
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Ditto of 



Ditto of 
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of «- 

2.3.4 


m 1 1 

~ 5 + 3m — 30m 3 




m 1 

— 2 ~ 2m" 




m 1 
— 3 _ 6m" 




+m 

11m 1 1 11m 1 — 10m 2 — 1 




30 3m 30m 3 — 30m 3 


(/) 
2. ...6 


(m 2 -l)(llm 2 +l) 

30m 3 
mil 1 

~ 7 ' 2m 6m 3 "*' 42m 5 




m 5 1 
2 4m 4m 3 




_ m _ JL . _L 

3m 6m 3 




15m 15 

+ 4 + 4^ 




4m 2 
+ _ 3" + 3m" 




— 6m 




191m 2 1,1 191m G -168m 4 -21m 2 — 2 




84 ' m ' 4m 8 ' 42m 5 84m 8 




(m 2 - l)(191m 4 + 23m 2 +2) 
84m 5 


. (A) 
2. ...8 


to 2 7 2 1 
~ 9 ' 3m 15m 3 ' 9m 8 30m 7 




m 7 7 1 
— 2 — 3m + 6m 3 — 3m 5 




_ 2m _7 _7_ J_ 
m m + 3m 3 3m 5 




28m 70 14 
+ ~3~ + 3m — 3m 3 




49m 49 49 
+ 5 + 3m - 30m 3 




_49m- — 
tn 




— 12m—- 
m 
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+72m 

2497m 24 98 4 I~~ 



90 m 30m 3 9m 5 30m 7 

_ 2497m 8 — 2160m 6 — 29 4m 4 -40m 2 -3 
~~ 90m 7 

_ (m 2 -l)(2497m 6 +337m 4 +43m 2 +3) 
~~ 90m 7 ' 

Heuce the formula which expresses the sum of the m— 1 inter- 
polated values, together with one-half of the two extreme values 

(m 2 -l)(191m 4 -|-23m 2 +2) 

60480m* ^ 

(m 2 — l)(2497m 6 +337m 4 +43m 2 +3) /7 , , N 
3628800m 7 W W 



Dividing by m we get the average of m values = 

dfflj (m 2 -l)(llm 2 +l) 
n 2 v ;_t " 720m" k 

(m 2 -l)(191m 4 +23m 2 +2), 



(V}_ 2-=l (5)+ («'-lXll""+l) ( . 
^ ' 12m 2 W+ 720m 4 W 



60480m 6 



(/) 



(m 2 — l)(2497m 6 + 337m 4 +43m 2 +3) 
+ 3628800m 8 ( ' )- 

When m is supposed to be indefinitely great, this becomes 

which is the same result as that determined by the integral calculus, 
and expresses the sum of an infinite number of ordinate spaces, or 
the curvilinear area contained between the ordinates V and V,, when 
the distance between them is unity. 

The formula (a) may be applied so as to include any number 

w-t-1 of primitive values V , V,, V 2 , V„. Denoting the 

respective coefficients by m, — B, D, — F, H, for the sake of brevity, 
the expression is 

2(V')o....,=m(V)-B(5)+D(rf)-F(/)+H(/0; 
and when extended, so as to comprise all the n intervals from V to 
V„, it gives 

2(V')o....„=mS(V)-B2(6)+DS(rf)-FS(/)+H2(A) 
=m2(V)-B{(a„)-(« )} + D{O_(r )} 

-F{(0-K)3+H{(y„)-(^o)3 • • • (»). 
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To employ this general formula, the values of (a ), (a n ), (b ), 
(b n ), &c., as ordinarily derived by a direct process of differencing, 
require a previous knowledge or calculation of certain values of V 
which precede V and follow V„. This, however, may be prac- 
tically dispensed with by supposing the differences to be extended 
upwards and downwards, so that the last order of differences shall 
retain the same constant value, or that the next following order of 
differences, and, of course, all succeeding orders, shall vanish or be 
disregarded. Thus, suppose the quantities V , Vj, V 2 , . . . . V„, 
to be differenced with a notation as follows : — 



v„ 
v, 
v 2 
v 3 
v 4 

V n _ 4 
V„- 3 

v„_ 2 

V„_i 

v. 


a 

«1 

a 3 

a'_ 3 
a _ 2 

a'-i 
a 


b 

h 
h 

v., 

b'-i 

b' 


c 
c 2 

C-2 

c'-i 
c 


d 

d, 

d'_2 

d'-i 
d' 


e 
e 

e 
b 


-l 


/ 
/l 

/'-l 
/' 


9 
9 


h 
h' 



Then we shall have 



«„+!=«'+ *'+ o'+ d'+ e'+/'+ /+ K 
.: (a„) =a'+$b'+ic'+W+y+W+y+Lh' 

a_! =a — b-\- c— d-\- e— J "+ g— A 

.-. (a )=a-$b+$c -&+&-$/+&-& 

c„_i=c'+ d'+ e'+ /'+ g'+ A' 
c n+l =c'+2d'+Se' + Af + bg' + 6A' 

••• (<v) =«/+K+2'+*/ + V+ S*' 

e_! =c-2d+3e— 4f+5g— 6A 

c +1 =e— rf+ e— /+ $r— A 

.-. (c ) = C -f rf+ 2e-|/+ 3 9~i h 



e n+ ,=e' + 3/' + 6y + 10A' 
••■(0=^ +*/ + #+ 7A ' 

e_! =e— 3/ + 6^— 10A 

e +l =e-2/+3g- 4A 

.■■(*)=«-*/ +fr- 7A 



<?„+i=/ + 4A' 
&)=• + **' 
^_i =g-ih 
9+i =ff—3h 
(g )=ff-lh. 
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And since 2(6) = («„) — (a ), 2(d) = (c„) — (c ), &c, these values 
give 

2(A) =(g-g)+ }{h' + A). 

The foregoing expression (5) for 2(V') ...» hence becomes = 

m2(V)— B(a'-a)— p(6'+6)-(£B-D)(c'-c) 
-(iB-|D)(rf'+rf)-(iB-2D+FXe'-e) 
_(lB-fD+#FX/'+/)-(iB-3D+|F-H)(/- i? ) 
-(p-fD+7F-^H)(A'+A); 

that is, after substituting the values of B, D, F, H, which are the 

coefficients of the formula (a), 

m 2 1 

2(V')o....„=m{(V 0+ V, .... + V„)-1(V„+V )}- 11 _ (a'-a) 

m2 ~W^ (^ 2 -l)(19^-l ) r , , (to 2 -1)(9to 2 -1) 



24m v y 720m 8 v ' 480m 3 

(»« 2 — l)(863m 4 -145m 2 +2) (m>—l)(275m*-61m* + 2) , ,, , ^ 

• (e — e) STTqoZS (/ +/ ) 



60480m 5 v ' 24192m 4 

(to 2 — l)(33953w 6 — 9247ot 4 +497w 2 -3) 



3628800m 7 
(w 2 -l)(8183m 6 -2617m 4 +197w 2 — 3) 



(9-9) 



1036800m* (A'+A) . . . (A).* 

By dividing by mn we get the mean value of V; and when m 
is made oo, the result multiplied by n gives the curvilinear area 
bounded by V and V„= 

(Vo-r-V, .... + V n )-i(V.+V )- l(a'-a)- l(6'+5) 

19 , • x 3 up i j\ 863 ., 275 

-72^ (c - c) -l60 (rf+ ' /) -60480 (e - fi) -24192 (/ + /) 

33953 , , 8183 tv x xs 



3628800 w Jy 1036800 

which is the formula usually employed for calculating quadratures 
where analytical integrations are impracticable. -See De Morgan's 

* This expression includes only one-half of the extreme values Vo, Y„. To obtain 
2 V"o....», including all themn + 1 values, primitive and interpolated, i(V„ + Vo) should 
be added, or the first term of the formula should be replaced by 

m(Vo + V, .... +V„) - *^ (V. + Vo). 
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Differential and Integral Calculus, p. 313, art. 181, where the series 
is given as far as (e'—e). 

In order that the investigation may be as complete as possible, 
we have carried it as far as the mean eighth order of differences, 
which is more ample than can be required, except for the purpose 
of ascertaining the magnitude of the terms rejected. 

For calculations in general, it will not be requisite to go beyond 
the fourth order, and, perhaps, in no case will it be desirable to 
proceed further than the sixth order of differences. Thus, dividing 
(A) by mn, we shall have for the average of mn values of (V), in 
which it has been considered sufficient to include the two leading 
terms of the several numerators, 

3(V')o....„ _ (Vq+V. • • • • +VJ-i(Vn+ Vo) 
mn n 

( rf_.) +1 (y+J)+g(rf-.)+i;(rf'+rf)+|S('-)+I^V+/) 



40 v 



5040 x 



2016 v 



12m 

The differences a, b, c, &c, may be found in terms of the 
primitive values by the simple algebraic operation of differencing, 
and the order of the symbols will be preserved by using the common 
rule of changing signs and adding, as in the following specimen : — 





Ai = « 


A 2 =4 


^3=0 


A 4 = d 


Vo 


-Vo + Vi 








V, 




Vo-2V l + V 2 








-V1 + V2 


V1-2V2 + V3 
\ r 2-2V3 + V 4 

V 3 -2V 4 + V 6 


-Vo + 3Vl-3V2 + V 3 




V'2 

v 3 

V, 


-V2 + V3 

-v 3 +v 4 


-Vi + 3V 2 -3V 3 + V 4 
-V2 + 3V 3 -3V 4 + V 6 


V -4Vi + 6V 2 -4V 3 + V4 

V 1 -4V2 + «V 3 -4V 1 + V 5 

&c. 




-v 4 +v 6 




&c. 


Vs 


&.C. 


&c. 







&c. 



On examining the process, it will at once be seen that the 
numeral coefficients, irrespective of sign, are made up according to 
the following table, in which the numbers in each column are 
deduced by respectively adding the number which stands in the 
preceding column to that which is immediately above it, until the 
said numbers arc severally exhausted. 
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1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


1 


1 


] 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 




1 


3 


6 


10 


IS 


21 


28 


36 


45 


55 


66 


78 


91 


105 


120 




1 


4 


10 


20 


35 


56 


84 


120 


165 


220 


286 


364 


455 


560 




1 


5 


15 


35 


70 


126 


210 


330 


495 


715 


1001 


1365 


1820 




1 


6 


21 


56 


126 


252 


462 


792 


1287 


2002 


3003 


4368 




1 


7 


28 


84 


210 


462 


924 


1716 


3003 


5005 


8008 




1 


8 


36 


120 


330 


792 


1716 


3432 


6435 


11440 




1 


9 


45 


165 


495 


1287 


3003 


6435 


12870 




1 


10 


55 


220 


715 


2002 


5005 


11440 




1 


11 


66 


286 


1001 


3003 


8008 




1 


12 


78 


364 


1365 


4368 




1 


13 


91 


455 


1820 






1 


14 


105 


560 








1 


15 


120 






























1 


16 
1 



These coefficients, which are here extended as far as the sixteenth 
order of differences for the purpose of reference, are evidently those 
given hy the binomial theorem, since by the operation of successive 
multiplication, these are obviously formed in precisely the same 
manner. 

Hence, according as n is odd or even, we have 



A»=- 



n— 1 




n—2 ■* 

^V 3 -&c. 



V.+&& 



(6).« 



Proceeding as far as the sixth order of differences, the following 
are therefore the several values : — 

a = -V +V, 
b= V -2V 1+ V 2 
c = _V +3V 1 -3V 2 +V 3 
<*= • V -4V,+6V 2 -4V 3 + V 4 
e = _Vo+5V,-10V 2 +10V 3 -5V 4 +V 5 
/= V -6V, + 15V 2 -20V 3 +15V 4 -6V 5 +V 6 
&c. &c. 

* If for Vo, Vi, Vj, &c, there be substituted their respective numeral coefficients, 
the formula (6) in either case will express the nth difference of those coefficients, which, 
being determined by a continuation of the process, must evidently be the central coeffi- 
cient of the 2nth order of the table. We thus arrive in a most simple manner at the 
curious summation 



i+ » !+ (»-V) 2+ 0- 



n-\ 
2 
= central coefficient of t^„ 

2n-l 2 »-2 
3 



=2»." 



re-2 



»+i 



y+&c. 



-r 1 (?)• 
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and these may, if required, be further extended by inspection of the 
preceding table. 

For certain special objects, such as the computation of areas, 
volumes, tonnage of ships, &c, it may be desirable to practically 
utilize the formula (C) by resolving a few particular cases, so as to 
express the results wholly in terms of V , V\, V 2 , &c. 

1. Let there be two intervals, or three given values ; theu n=2, 
and 

a=-V +V„ a'=-V,+V 2 , 6=6'=rV -2V, + V 2 

.:a'-a= V _2V, + V 2 
6' + 6=2(V -2V, + V 2 ) 

and the average or 2m values = 7 

+ iaS» — • • • (Ca) - 

2. Let there be three intervals, or four given values; then 
n=3, and 

a = - V + V, b = V -2V, + V a 

a'=-V 2 +V 3 i'=V,-2V 2 + V 3 

c=c = -V +3V,-3V 2 + V 3 

.■.a'-a=V -V,-V 2 + V 3 
i' + i=V -V l -V 2 +V 3 
c'-c=0 

and the average 01 6m values = — - 

o 

, (V +V 3 )-(V, + V 2 ) 



24m 2 
3. With four intervals, or five values, n=4; 

a'-a= V - V, • _ V 3 + V 4 
b' + b= V -2V l + 2V 2 -2V 3 + V 4 
c'-c= V -4V, + 6V 2 -4V 3 +V 4 
rf' + rf=2(V -4V I + 6V 2 -4V 3 + V 4 ) 

and the average of 4m values = 

7(Vo + V 4 ) + 32(V, + V,) + 12V, 
90 
, 7(V 0+ V 4 )-16(V, + V 3 ) + 18V 2 



(-3). 



144m 2 

4. With five intervals, or 6ix values, n=5 ; 

a--a=V - V, • • - V 4 + V„ 

i' + 6=V -2V,+ V 2 + V 3 -2V 4 + V 5 
e'-c=V -3V 1 + 2V 2 + 2V 3 -3V 4 + V 5 
rf' + fi ?=V -3V 1 + 2V 2 + 2V 3 -3V 4 + V 5 
e— e=0 



(««)• 
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and the mean of 5m values = 

19(V +V 5 ) + 75(V, + V 4 ) + 50(V,+ Va) 
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+ 



288 - 

5(Vo + V»)-9(Vi + V4) + *(V g + V a) 
144m 2 



(<■.)• 



With six intervals, or seven given values, n=6; 



- a = V _ V, 



-2V,+ 



» 



* + 



* - v,+v 8 

V 4 -2V 5 + V 6 



b'+b: 

c'-c= V -3V,+ 3V 2 - 2V 3 + 3V 4 -3V 5 + V 6 
d' + d= V -4V,+ 7V 2 - 8V 3+ 7V„-4V 5 + V 6 
e' _,= V -6V 1 + 15V 2 -20V 3 + 15V 4 -6V 5 + V (! 
/' + /=2(V -6V 1 + ]5V J -20V 3 +15V 4 -GV 5 +V C ) 

and the mean of 6m values = 

41(Vq+V6)+216(V,+V 5 ) + 27(V 2 +V4)+272V 3 

840 
1 57(V + V 6 ) - 432( V, + V,) + 675( V 2 + V«) - 800V 3 



+ 



4320m 2 
/ 



37/ 



, , for the purpose of 



+ 



• • («•«)•• 



or (with the small variation g^ — ^q<>o ot 2 > 
simplification) = 

(Vo + V6)+5(V. + V») + (Va+V4) + CV3 
20 

4(V o+V 6 )-7(V 1 + V.,) + 4(V 2 + V 4 )-2V 3 
144m 2 

It may be of some use, in forming practical notions, to briefly 
advert to the geometrical signification of the set of formulae just 
given. In the diagram, as before, 
the given series of quantities V , 
Vi, V 2 , &c, are represented by 
the darker ordinates, the upper 
extremities of which are con- 
sidered to be sufficient to deter- 
mine the contour of the curve 
line P P, P 2 &c. If these 
extremities be also connected by 
straight lines, PP,, P,P 2 , P 2 P 3 , 
&c, a polygonal figure will be 
formed. Suppose that m = 5, » = 6, or that four ordinates are 

* An expression equivalent to the first part of this last formula, for finding the area 
of a figure, was given by the late Thomas Weddle, in the Cambridge and Dublin Mathe- 
matical Journal, New Series, vol. ix., p. 79. 




Mm Mi M» MJ M4 MS 
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to be interpolated in each of six intervals, as shown in the diagram. 
Then, if the small subdivisions (Mm) of the line of abscissa? be 
taken as unity, we shall have MM]=m, and the total abscissa or 
base line MN=m». Also, if we consider the polygonal figure 
bounded by the extreme ordinates and the base, we shall have 

Polygonal area=m{(V + V 1 . . . . + V n )-l(V 0+ V n )}. 

Consequently, the corresponding average of the given quantities 
Vo> "Vi, . . . . V n , estimated as before described, is 

(V„-fV,....+V,,)-l(V +V„) = Polygonal area 
n Base 

Again, if after the insertion of the interpolated ordinates the 
upper extremities of the entire set be consecutively joined, an inter- 
polated polygonal figure will be formed, the area of which will be 
expressed by the sum of all the ordinates after deducting one-half 
of the extreme values, that is, by S(V') ... n ; 

S(V') „ Interpolated polygonal area 
mn Base 

It is hence evident that, if the formula (C) be multiplied by the 
base MN (=mn), the respective terms of the product will be thus 
geometrically represented : — 

Interpolated polygon = Primitive polygon 

— Curvilinear segments 
+ Interpolated segments. 

Also, if the particular formulae (c 2 ), . . . (c 6 ), be severally mul- 
tiplied by the base MN(=mn), the geometrical signification of the 
terms will become 

Interpolated polygon=Curvilinear area 

+ Interpolated segments. 

In these relations, it will not pass unnoticed that the several 
segments must be considered as positive or negative, according as 
the convexity of the curve may happen to be presented downwards 
or upwards. 

It may be further observed, that the term which expresses the 
segments contains the reciprocal of m? as a factor ; from which we 
draw the practical conclusion, that if a number of small segments 
of nearly equal length be cut off from a given curvilinear area, so 
as to reduce it to a polygon, the value of the sum of these segments 
will vary inversely as the square of their number. 

By way of showing the degree of accuracy of the formulae (c 5 ), 
(c 6 ), we shall employ them in finding the quadrature of the circle. 
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Let CD be an ordinate bisecting the radius OA of the quadrant 
AOB, and suppose four additional equidistant 
ordinates to be interposed between BO and CD. 
To obviate fractions, take 0A=0B=10; then 
OD = 5, and the six abscissae, estimated from 0, 
will be respectively 0, 1, 2, 3, 4, 5. The cor- 
responding ordinates are, therefore, V =10, 
V, = 1/100-1=31/11, V 2 = 1/100- 4=41/67 V 3 =V / T00 
= 1/91, V 4 = 1/100- 16=21/21, and V s = 1/100-25= 5 ^W. 
Hence, by the first part of the formula (c 5 ), the mean ordinate of 
the curvilinear area BD is found to be 

(V)= JL { 1 9(2 + VW) + 1 5(3 V lT+ 2 \/21) + 10(4 V 6 + V9l) ] . 

Also the area B0DC=0D x (V), and the area of the triangle 
COD =0D x £CD =0D x £V 5 ; .-. the sector B0C=0D { (V) -$V 6 } . 
Now, as BC is an arc of 30°, the area of the quadrant AOB is three 
times this sector, and if we divide by r 2 =100, to reduce the scale 
of measurement to radius = unity, the result is 

_L Jr\n_i.v 1 - 38 + 5(2\/9i + 8\/6"+6^2l + 9\/Tl-25v^3) 
20 1 1 } 2 Vs J- 384 

38 + 5( ^364 + VMi+VTb6 + V89T- a/187 5) 

- 6.8.8 ~~ 

= -7853973, 

which is very easily found from a table of square roots. This 
computation diners from the true value ('7853982 &c.) by only an 
unit in the sixth place of decimals, and gives a nearer approximation 
to the circle than would be obtained by taking an inscribed polygon 
of 2,000 sides. 

To determine the same area from seven ordinates, or six intervals, 
it will be requisite to interpose five ordinates between BO and CD ; 
and fractions will be avoided by first taking 0A=0B = 12. Thus, 
OD = 6, and the seven ordinates are respectively 



Vo = 









12-0000000 


V, = 


V1U- 


-1 = ^143 


11-9582607 


v 2 = 


\/144- 


-4 = \/140 


= 


11-8321596 


v 3 = 


VlU- 


-9 = VUb 


= 


11-6189500 


v 4 = 


V144- 


-16= v'm 


= 


11-3137085 


v, = 


V\U- 


-25= V119 


= 


10-9087121 


v 6 = 


^144- 


-36= \^108 

••• v +y 6 


== 


10-3923048 




22-3923048 






v,+v 5 


= 


22-8669728 
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and by substitution in the first part of the formula (c 6 ) the mean 
ordinate is 

(V)=ll-4793368. 

Also, as before, the sector B0C = 0D |(V)— ^V 6 |, which, mul- 
tiplied by 3 and divided by r 2 =144, the area of a quadrant, to 
radius=unity, comes out '7853981, which is only an unit in error 
in the seventh place of decimals, and no part of the calculation has 
been taken beyond that place.* 

Having completed the theory of summation, as originally pro- 
posed, and endeavoured to elucidate descriptively and practically 
the principles upon which it rests, we shall now take up the subject 
in a different manner, and show how the general expression (5), 
p. 307, may be transformed into a new formula, which, after some 
slight modification, is applicable, in many cases, to a yet more simple 
and easy calculation of areas or volumes, and is also most especially 
adapted to expedite by a brief and elementary process the com- 
putation of the values of annuities that involve the contingencies of 
any number of lives. 

Restoring the values of the coefficients in terms of m, the equa- 
tion (5) gives 






(m 2 -l)(191m'< + 23m 2 + 2 i. . , ,, 
+ 60480^ }(*o)-(0} 



&c. 

=S(V) 

f»» a -l, , llm 4 -10m 2 -l. . 191m 6 - W8m 4 -21m?-2, s ] 

+ i-12^- (0o) 7M (r ° )+ 60480^ ("o)-&c.f 

m*-\, . 11m 4 - 10m 2 - 1. 191m 6 - 168m 4 -21m 2 -2 . . „ 1 
TM M 720^ (C " )+ 60480^ 0O-&cJ. 



"I 



Now, if we suppose a value Y x to be taken near to V , and 
divide the increment V x — V by x, the quotient, after x is made 



* For the ratio of the circumference of a circle to its diameter I have, by means of 
an empirical process, constructed the concise numerical expression 

which is extremely simple considering the degree of accuracy attained. It is correct to 
no less than ten decimal places of figures. 
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dV 
=0, will express the value of the differential ratio -*- at that point. 

Hence, from the expression (a) in terms of the central differences 
(article on "Interpolation," p. 68) we thus find, at the initial 
point, 

Similarly at the terminal point, where x=n, we shall have 

dV\ , . 1 . . , 1.4 , , „ 

Therefore, by substitution, we get 

m y Jnr 12m? \\dx J \dx J J 

m*-l . 23m«-21m?-2 

~ 720^ i^-W} + 60480m' ~ {W" W}-&c, 

where the latter terms, being respectively less than =^ { (c ) — (c„) \ , 
noon { ( e o) — { e n) }> &&> may be rejected as insignificant. 



( 



The expression thus reduces down to the brief and neat form 

'■Jo \ <*»/».' 



w .A^-r 12m2 ^^ ; ' 



JM 2 — 1 

=S(V)+ — -5- (tan e - tan 0„) 



(A 



where 8 denotes the angle which a tangent to the curve makes with 
the line of abscissae. As the interval between the interpolated 

ordinates is — , this expresses the area of the interpolated polygon, 

and the first term S(V) expresses that of the primitive polygon. 

* By successively differentiating the expression (a) referred to, and employing the 
initial values of the odd differential coefficients, that is, the values when a=0, the 
formula may be completely resolved into the following remarkable development : — 



S(VQo....„ _ 



m 



„.„. m 2 -l ( sdV\ sdV\ ) 

n*-\ \ f<PV\ _ (-<PV\ 
120m* ( V<W \<h?)„ 



m*- \ 
~ 720m 

ro 6 -l I /<PV\ _ s<PV\ ) 
+ 30240m« I \d^J ydat'n) 
- &c. &c. 
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Also, by making the number m indefinitely great, the curvilinear 
area is 



^m.-mn 



= 2(V)+— (tane -tane„) 



(ay 



And from these expressions we further conclude, that the sum of the 
curvilinear segments = yn (tan0 o — tan#„), and the sum of the 

interpolated segments = fn~i ( tan ^o — tan 9„) , it being remem- 
bered that tan is positive when the curve ascends, and negative 
when it descends in a forward direction. 

The formulae (d) and (oY) are remarkably simple, and suggest 
the following geometrical statement : — 

In the diagram P, Q, are the initial and terminal points of the 
curve; V , V], are the first two, and 
V„_i, V„, the last two of the given 
ordinates. Let tangents PT, QT", be 
drawn touching the curve at the two 
extreme points, P, Q, meeting the 
adjacent ordinates, V,, V n _ u pro- 
duced if necessary, in T, T'; and let 
Pm, Qw', be drawn parallel to the 
horizontal line of abscissae, meeting the 

same ordinates in m, m'. Then (P) denoting the polygonal area 
MP?! .... QiQN, and the intervening distance MM, being not 
necessarily considered as the linear unit, we shall have 

Polygonal area (P)=2(V).MM, 

= {(V +V,. . . . +V„)-i(V 0+ V„)}.MM 1 . 

After m— 1 ordinates are interposed between each of the n 

intervals, 

MM 
Interpolated polygon (F)=2(V).- 



P 1 


A 


T T 

m 

r 

m 


r^ 


\ 


\ 




%-! 



v„. 



M M, 



3T, S 



=(P) + 



=(P)" 



m 

m 2 -l 

12m 2 

m 2 -l 

6m 2 



(Tm-t-TW).MM, 
(T + T'). 



* Here the interval or distance between the given ordinates is supposed to be unity. 
If the distance be o the curvilinear area will be 

o 2 
=a.2(V)+ p^ (tan 9o - tan 0„)- 
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And, after the interpolation of an indefinitely great number of 
ordinates, 

Curvilinear area=(P)+ — (Tm+T'w^.MM, 

1 £i 



=(P)+g(T+T'), 

where T, T', denote the respective areas of the tangential triangles, 
TPm, T'Qm'. 

Hence, the sum of the areas of the curvilinear segments cut off 

T + T 

by the sides of the primitive polygon = — ^ — , 



and the sum of 



T + T' 
those cut off by the sides of the interpolated polygon, = fi 8 . 

The lines Tm T'm', and the tangential triangles T,T', must, of course, 
be accounted as negative in value whenever the points T, T', may 
happen to fall below instead of above the horizontal lines which 
pass through the extreme points P, Q. They are, in fact, positive 
or negative, according as the tangents ascend or descend from the 
point of contact. 

From the generality of the investigation by which these results 
have been obtained, they are evidently independent of the existence 
of any law of sinuosity or flexure. The only essential requisitions 
are, continuity of outline, a nonvertical position at the extreme 
points, and a sufficiently narrow interval between the ordinates to 
include peculiarities. That is, with the knowledge of the initial 
and terminal directions, PT, QT', the intermediate points, P u P 2 , 
&c, should be sufficiently numerous to indicate the track of the 
curve. 

A somewhat curious conclusion may be drawn, if we conceive 
the case of a curve having contrary flexures, which bring its course 
back again to its primitive direction, 
as, for example, in the annexed 
diagram. The tangents, PT, QT, , 
are here parallel, and the tangential 
triangles are equal in magnitude 
but contrary in algebraic sign, so 
that their sum is zero. Hence, in 
this particular case, the polygonal 
and curvilinear areas are equal. 
This result is, in some \ measure, 

explained by the apparentjact^that, according as the convexity of 

z 2 
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the curve is presented upwards or downwards, the curvilinear seg- 
ments are severally added to or cut off from the polygonal area; 
and it is remarkable that these alternations, when taken in com- 
bination, should become practically balanced. 

Independently of the mathematical investigation, a little con- 
sideration will ,elicit a popular explanation of the philosophical 
nature of these geometrical deductions. Thus, when the curve has 
but little curvature, and alters its course gradually, the segments 
cut off by the chords are severally smaller in magnitude, but as the 
curve is longer in effecting a given change in its direction, a greater 
number of segments will necessarily intervene. On the other hand, 
when the curve has more curvature, and turns more suddenly, it 
will be perceived that the segments are comparatively greater in 
magnitude, but, as the curve sooner acquires its new direction, they 
are fewer in number. In the two cases, therefore, when the values 
are estimated between given tangential directions, the sum of the 
segments may be regarded as having practically the same value. 

In order to adapt the formula to the computation of life annui- 
ties, let the unit of interval be taken between the interpolated 
instead of the primitive values, which will now be, V , V m , 
Y 2m , .... V„. Also as an interval of time (t) expressed in years 
or parts of a year, is measured accordingly, the symbol x, denoting 

the abscissa, should obviously be replaced by — . Thus, the formula 

(d) becomes 

But when the annuity is completed, and the series representing 
it is entirely run out, the final number will extend beyond the 
limiting or highest age of the table of mortality, and we shall then 

evidently have V„=0, ( -5- J =0; 

, 2( v,™; ( v )+ ^(g). 

Let now p denote the yearly interest on £1 ; 

v the present value of £1 certain receivable a year hence ; 
l the number living at the given age according to the 

table of mortality; 
I or l t the number living after t years have elapsed; 
V or V* the present value of the fth annual payment of 

£1 receivable by the annuitant on the expiration 

of t years; 
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then \= T v' 

'o 

•log V=log l-\-t log v— log / ; 
and, differentiating with respect to t, 

\ d\ 1 dl , 

— • —r- = r • — + log V 

V dt I dt 6 

= _p-log(l+p) 

= -(r« + 3), 
where ju denotes the force of mortality at the given age, and 8 the 
rate of interest when it is assumed to be payable momently. The 
latter quantity (§), as it plays an analagous and equally important 
part in the effective diminution of the value of V, we shall, with 
like propriety, designate the force of discount. 

Now, if A denote the required value of the annuity, and we 
calculate 

S=V m +V 2ra + V 3ra + &c, 

we shall evidently have 

S(V')=iV +(V l+ V 2 + V3 + &c.) 
= i -f-A 
(V) =Iv +(v„+ V 2ra+ V 3 „+&c.) 

By substitution, the formula hence becomes 

i+A=^+s)-^ r -( A «+a) 

...A=«8+==- 1 -I^=i0'+*) • • • P), 

and it will be perceived that S is the present value of an annuity 
payable only after intervals of m years, which, when m is taken at 
a wide interval, will consist of very few terms, and involve com- 
paratively but little calculation. 

It may be observed as the peculiar and valuable feature of the 
formula (D), that, irrespective of the number m, the data 

correctly fix the initial position and direction of the annuity curve. 
Thus, when we have to consider the assumption of the number m, 
in connexion with the important adjustment alluded to, it will 
appear that the subsistence of very few of the distributed ordinates, 

* In mathematical analysis hyp. log. is always understood, the common logarithm 
being exclusively employed as a subsidiary in numerical calculations. 
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V m , V 2m , V 3m , &c, will be sufficient to impart to the curve its 
specific form and character, and thereby ensure a requisite degree 
of accuracy in the result of the computation. 

The formula will apply equally to a joint annuity on any number 
of lives, if only fx + fi + &c., be substituted in place of fi, according 
to the following theorem : — 

Theorem. — The force of mortality of a joint existence is equal 
to the sum of the forces of mortality of the constituent lives. 

The demonstration may be readily inferred from the differential 
relation 



d(W &c.) 

//' &c. 



dl dl' „ 



MORTALITY 
ONLY 



M, 



DISCOUNT 
ONLY 



M 



T 



MORTALITY 
AND 

DISCOUNT 



The precise adaptation of this theorem to the method before us 
leaves nothing to be desired, as it at once brings the calculation of 
a joint annuity on any number of lives under the same elementary 
process as that for a single life. We may here briefly state the 
corresponding geometrical relations. 

If, in the first diagram, the ordinates PM, pM u &c, be sup- 
posed to represent the pro- 
portion of survivors after 
0, 1, &c., years, the curve 
Vp, &c., will be the curve 
of mortality; and if Vm be 
drawn parallel to MM, and 
Vt a tangent to the curve 
at P, meeting the second 
ordinate ^?Mi produced in 
m and t, then^will mt repre- 
sent fi the force of mortality at P ; and if the curve of mortality be 
that of several lives jointly, it will, by the theorem, be the sum of 
the individual values of fi. 

Again, in the second diagram, if the ordinates, PM, /M„ &c, 
be supposed to represent 1, v, v 2 , &c, the present values of the 
payments of an annuity certain, then the logarithmic curve Pp', 
&c, will be that of discount, and it will differ from the former by 
its convexity being uniformly presented downwards. Drawing the 
lines Pm', Vtf, as before, m't' will represent the value of 8 the force 
of discount, which is invariable and independent of age. 

Lastly, in the third diagram, if the ordinates represent the 
values of the proportion of the survivors discounted for present 
values, the curve PP! , &c, will be that of the annuity payments at 



M, 



M 



M, 
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the proposed rate of interest, and mT the total force of decrement 
will, according to the theorem, be the sum of those of mortality and 
discount, shown in the two former diagrams. 

In the several diagrams, only one yearly interval is exhibited, 
the others not being required for our present object. It will be 
observed, that in each case the force of mortality, discount, or of 
total decrement is the dip or vertical descent of the initial tangent 
when it has reached the next yearly ordinate. 

Resuming now the general case of a joint annuity, if p t , p' t , &c, 
be the probabilities of the respective lives surviving t years, ptp't 
&c, or ( pp' &c.)( will express the probability of their joint existence 
after / years. Thus we shall have 

v » =(PP' &c.) m .« ra 

Y 2m=(PP &C.) 2m .« 2 " 

&c. &c. 

A=m(V m+ V 2ra +&c.)+ ?y -^^(n + fi' + &c + 2) . . (D'). 

It is another material advantage in the present method that the 
formula may be employed with equal facility and without any addi- 
tional labour, when the joint annuity is held on a number of lives 
which, taken separately, are classed under different tables of mor- 
tality. We have merely to take the probabilities of survivorship 
(p, p', &c.) and the forces of mortality (ju, //, &c.) from the respec- 
tive tables, and the form and extent of the calculation is precisely 
the same. 

The interval m may, in most cases, be extended up to about 
20 years. In general, however, the best intervals to adopt will be 
5, 7, or 11 years. For these intervals, the coefficients, &c, are 
convenient integers, and the special formulae are 

m= 5 years, A= 5S + 2— 2(/x + S) 
»i= 7 „ A= 7S + 3- 40 + S) 
w=ll „ A=llS + 5-10(/u+g). 

To practically elucidate the method, tables are here given 
showing the values of the force of mortality (/u) at each age, accord- 
ing to the Experience, Davies' Equitable, and Milne's Carlisle 
tables, and the force of discount (8) at several rates of interest. 

The values of the force of discount have been calculated from 
the formula S = hyp. log (1 + p) . It is nearly equal to \ (1 — v + p), 
that is, the average between the discount and interest of £1 for 
one year, and it is in reality that rate of interest which, if accumu- 
lated momently, would be equivalent to that which actually exists. 

The force of mortality has been found by dividing half the sum 



324 



On Interpolation, Summation, fyc. 



[JULT 



of the decrements immediately preceding and following by the num- 
ber living at the given age. The formula (8) might have been made 
use of as the correct numerator, but with small numbers no increase 
of accuracy would be gained. Any irregularities of progression 
must be attributed to the imperfect adjustment of the original tables. 
Force of Mortality (fi). 



Age. 


Experience. 1 


[ilne's ] 
arlisle. Eq 


)avies* 
uitable. 


Age. 


Experience. 


Milne's 

Carlisle. 


Davies' 
Equitable. 


10 


•00678 


0048 


0039 


51 


•01655 


0139 


•0189 


11 


•00680 


0047 


0039 


52 


•01757 


•0149 


•0204 


12 


•00682 


0049 


0041 


53 


•01868 


•0158 


•0219 


13 


•00685 


0051 


0043 


54 


•01989 


•0167 


•0230 


14 


•00690 


0054 


0045 


55 


•02120 


•0176 


•0238 


15 


•00695 


0059 


0048 


56 


•02263 


•0186 


•0250 


16 


•00700 


0065 


0052 


57 


•02418 


•0201 


•0262 


17 


•00706 


0068 


0056 


58 


•02585 


•0228 


•0276 


18 


•00712 


0070 


0060 


59 


•02767 


•0265 


•0290 


19 


00719 


0070 


0064 


60 


•02970 


•0313 


•0302 


20 


•00727 


0071 


0066 


61 


•03195 


•0352 


•0311 


21 


•00736 


0070 


0067 


62 


•03442 


•0373 


•0325 


22 


•00745 


0070 


0069 


63 


•03712 


•0386 


•0343 


23 


•00754 


0070 


0072 


64 


•04008 


•0398 


•0363 


24 


•00764 


0071 


0074 


65 


•04332 


0413 


•0384 


25 


•00779 


0072 


0077 


66 


•04686 


•0427 


•0408 


26 


•00786 


0074 


0079 


67 


•05073 


•0444 


•0434 


27 


•00798 


0076 


0084 


68 


•0.5495 


•0464 


•0463 


28 


•00810 


0083 


0088 


69 


•05950 


•0489 


•0494 


29 


•00824 


0093 


0093 


70 


•06443 


•0516 


•0525 


30 


•00838 


0100 


0096 


71 


•06980 


•0567 


•0560 


31 


•00854 


0102 


0099 


72 


•07563 


•0653 


•0598 


32 


•00870 


0102 


0104 


73 


•08195 


•0756 


•0642 


33 


•00887 


0101 


0107 


74 


•08882 


•0874 


•0693 


34 


•00905 


0102 


0108 


75 


•09631 


•0973 


•0738 


35 


•00923 


0103 


0112 


76 


•10442 


•1043 


•0789 


36 


•00943 


0105 


0115 


77 


•11326 


•1111 


•0849 


37 


•00963 


0108 


0118 


78 


•12295 


•1146 


•0918 


38 


•00984 


0111 


0122 


79 


•13350 


•1203 


•1011 


39 


•01007 


0116 


0124 


80 


•14496 


•1280 


•1125 


40 


•01030 


0125 


0125 


81 


•15739 


•1362 


•1256 


41 


•01054 


0135 


0127 


82 


•17083 


•1476 


•1408 


42 


•01081 


0142 


0128 


83 


•18546 


•1573 


•1589 


43 


•01113 


0146 


0132 


84 


•20157 


•1682 


•1790 


44 


•01154 


0148 


0139 


85 


•21958 


•1820 


•2009 


45 


•01202 


0149 


0143 


86 


•24025 


•2030 


•2238 


46 


•01260 


0149 


0145 


87 


•26419 


•2280 


•2448 


47 


•01326 


0148 


0150 


88 


•29247 


•2478 


•2613 


48 


•01398 


0144 


0157 


89 


•32672 


•2486 


•2706 


49 


•01476 


0139 


0165 


90 


•36847 


•2676 


•2769 


50 


•01562 


0136 


0176 











Force of Discount (2). 





f 1 


•00995 


Interest 


H 


•01489 


per Cent. 


2 


•01980 




in 


•02469 



Interest 
per Cent. 



I 3 
34 

4 
I 44 



•02956 
•03440 
03922 
•04402 



Interest 
per Cent. 



•04879 
•05827 
•06766 
•07696 
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Example 1. — Required tbe present value of an annuity payable 
during the existence of a life aged 40, according to the experience 
table of mortality and interest at the rate of 3^ per cent. 

Take m=7 years, and from Y t —p t .vf we find 





v 7 


= -7253 




V 14 


= -5089 




v 21 


= -3351 




v 28 


= -1959 




V35 


= 919 




V42 


= 291 




V49 


= 44 




v 56 
s.. 


= 1 




..1-8907 






7 




13-2349 


H -01030 




+ 3 


a 3440 












16-2349 




•04470 x 


-4= 


-•1788 



A.... 160561 
Otherwise, take m = ll years, then 

V„ = -5957 
V 22 = -3132 
V 33 = -1177 

V 44 = 187 
V 55 = 2 



S.... 1-0455 
11 



11-5005 
fi -01030 +5 

8 3440 

16-5005 

•04470 x -10= — -4470 



A.... 16-0535 

Jenkin Jones, in his Experience Tables, entitled New Rate of 
Mortality, gives the value 16055. 

Example 2. — Required the value of a joint annuity on three 
lives, aged 35, 60, 65, according to the Carlisle table, at 5 per cent. 



326 On Interpolation, Summation, tyc. [July 

Take m=5, then from Y t =(pp , p // ) t .v t we find* 







v. 


= -4888 






v,„ 


= -1980 






v lg 


= 573 






V, 


= 110 






v M 


= 12 






V30 


= 1 




•0103 


V35 

s. 


= 


M 


...•7564 


/*', 


313 
413 




5 


/*" 










3-7820 




•0829 




+ 2 


s 


, 488 






c 1 . 






— — 




5-7820 




1317 


X-2: 

A. 


= -•2634 




...5-5186 



Milne, in his Treatise on Annuities and Assurances, Table L., 
p. 676, gives the value 5'519. 

By means of the formula (d), we can deduce the correction that 
should be made in the value of a life annuity when each yearly 
payment is made in m instalments. Let V l3 V 2 , V 3 , &c, be the 
present values of £1 receivable by the life, respectively after 1, 
2, 3, &c., years, so that 

A=V 1 +V, + V, + &c, 
and the value of the annuity, when each £1 is payable in m instal- 
ments of — each, will be 
in 

A'= — (V' 1 + V' 2 + V' 3 + &c.) 
m 

Also, since V =l, we shall evidently have S(V)=^ + A and 

dV 
2(V') = J + biA'; and, beyond the limiting age of the table, -7- =0. 

Therefore (d) gives 

i + ™ A ' _i 1 a 1 r^-\ (dV\ 
m _2+ + 12m* \dt) 

* The values of logp< are easily got out from a table of log I, by placing a card over 
log? and another under log / ( , so as to obtain their difference without transcribing them; 
and if for any one of the lives a table of log D is available in place of log /, the factor v' 
will be thereby included, which will so far abbreviate the calculation. The logarithms 
taken to five places will be quite sufficient. 
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,, , m—l m?—\ , ., ... 

•■•A=A + — - 1S ^0. + «). • -(9) 

where the positive term is the polygonal correction, and the nega- 
tive term the correction for curvature, the latter of which has 
generally been disregarded. 

Hence, when a life annuity is payable 

Half-yearly, A'=A+ \ - — (/i + 3)' 
5 



Quarterly, A'=A+ f- — (ji + S) 



1 



(10). 



Momently, A = A + £ - — (fi + 8) I 



The present value of a life premium of £1 payable yearly, being 
that of an annuity payable in advance, = 1 + A. Add an advance 

payment = — to (9) and the present value of a life premium of 

£1 per annum, payable in m instalments, is 

m + 1 to 2 — 1 , 

ni — 1 »» 2 — 1 , „, 

Therefore if tt denote the ordinary annual premium for a life 
assurance, the same when payable in m instalments will be 

t(1 + A) 



.v m—l m' 2 — 1 , . 

(1+A)-— — — (ji + i) 

2m \lm l 



sfm — 1 mi — 1^ ^ 



(ii). 



According to the formulae hitherto given, the value is somewhat 
underrated, more especially at the higher ages.* 

The formula (a") will also serve to determine, from a given 
table of mortality, the mean duration of life after having attained a 
given age. This is more commonly called the " expectation of life," 

* In official practice it is a safe and convenient rule to take the yearly premium a year 
higher when it is to be paid quarterly, and half a year higher when it is to be paid half- 
yearly. This will rather more than cover the true calculation, as it should do to provide 
for the cost of additional receipts, &c. 



=H-^V^ -• • • -(12). 
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and, in the tables usually given in works on life contingencies, it is 
calculated on the supposition that the annual decrements are 
equally distributed throughout each year. As this assumption is, 
in fact, equivalent to the substitution of a polygon in place of the 
true curve of mortality, it must evidently involve some error, and in 
summing the terms this must accumulate from year to year. If, 
in the mortality table, l , l u l 2 , Sec, denote the successive numbers 
living, commencing with the given age, and these be taken as 
ordinates and the number of years elapsed as the corresponding 
abscissae of the curve, then the formula (d') will express {fldt) the 
total number of years of life eventually passed through by the 
persons originally represented by the number l . Hence, observing 

that at the end of the table -j —0, and dividing by l , we find the 

mean duration of life to be 

_S(Z) 1_ }_(dl\ 

-~T; + 12 'T [di) 

1 , li + h + h+&C. 

k 

It is, in fact, the same as the value of an annuity payable momently, 
and supposing money to have no interest. Thus it appears that 
the tabular expectations of life, usually given, are in excess of the 
truth by one-twelfth of the force of mortality at the respective 
ages. The values of this force, of which a table has been given, 
will show that the requisite correction is decidedly appreciable at 
the advanced ages, but the error is not of much consequence, 
as the mean duration of life is an element that is only made 
use of in cursory comparisons, and is not required in professional 
calculations. 

Reverting, for an instant, to the application of the general 
subject of investigation, we have found it to resolve itself into the 
determination of a convenient practical relation expressing with 
sufficient accuracy the value of the segments which make up the 
difference between a curvilinear area and that of an inscribed 
polygonal figure, the points of which are marked out by a given 
series of equidistant ordinates. Before closing this article, it may 
not be uninteresting to subjoin a rather neat geometrical represen- 
tation of the corresponding relation when the several sides of the 
inscribed polygon are equal in length. 

Let PqP^j, &c, be a polygon of equal sides inscribed in a 
curve, and, taking P as origin, let the length of the arc of the 
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curve, to any point p, =s; the angle which the tangent to the 
curve makes with the initial tan- 
gent = 0; the radius of curva- 
ture, in the vicinity of the point 
p, =R; and the chord P Pi> or 
side of the polygon, =c. Then 
by disregarding, as inappreci- 
able, certain small terms of the 
fifth and higher orders, it may be 
shown, by the differential calcu- 
lus, that the curvilinear segment 
between P and P! may be estimated as equivalent to that of a 
circle having the average curvature between those points. Now, if 
ds, d9, denote corresponding small increments at the point p, we 

ds 
shall evidently have ~RdQ=ds, or R= -m; and therefore the average 

value of R, for all points of the curve lying between P and Pi, is 

where &, $0, are the values estimated between the extreme points 
of the segment. 

Hence, treating the segment of the curve from P to Pi as 
equivalent to that of a circle whose centre is and radius = (R), 
and observing that the angle P o OPi = S0, we shall have 
c 2 =(R) 2 ch =2(R) 2 (l-cosS0) 

Also, sector P o 0P 1 =|(R)^=KR) 2 ^9 



-»M 



&c. 



triangle P 0Pi=KR) 2 sin id 
segment P Pi=i(R) s (3e— sin 59) 



1 2. 
30* 



~12 V 4 - 5 ' 4 - 5 - 6 -7 



-&c.) 

+ &C.V . 



.7 
— &c. 



(13) 
(14) 
(15). 



* A farther extension of this formula will show a remarkable approximation of the 
terras to a geometrical series. 
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In this last expression the reciprocal fraction within the paren- 
thesis is so small that in the process of summation it may be con- 
sidered as a constant depending on an average value of $9, which 
again depends on the total number of sides of the polygon. Thus, 

if n denote the number of sides, the average value of $0= — 

n 

Q(") 30 

(where tt=3-141592, &c); and .-. ~m = —.^=-76^; and (3) 
hence becomes 



Segment=^[l+. 7Gw2 



M < 16 > 



Therefore, summing for any number of consecutive segments, 
we get 

(2 segments) .... m = ~G m \^l + . 76w> _ 1 ] .... (17), 

which is easy of calculation, and suggests the following compendious 
geometrical relations, which will hold generally, whatever changes 
may occur in the flexure of the curve. 

With centre P and radius equal to the side of the polygon 
describe a circle, and draw the tangential radius P T , which will 
indicate the direct course of the curve at the point P . Similarly, 
if we conceive tangents to be drawn to the curve at the remaining 
points, P], P 2 , &c, and radii P T], P T 2 , &c, to be drawn parallel 
to these tangents, the complete set of radii thus determined will 
respectively indicate the direct course pursued by the curve at the 
corresponding angles of the polygon. 

Then, the respective areas of the curvilinear segments, P Pj, 
PiP 2 , &c., will be practically represented by one-sixth of the circular 
sectors contained between the successive radii P T , P Ti, P T 2 , &c. 
Also, it follows that the sum of any number of consecutive curvi- 
linear segments will be represented by one-sixth of the circular 
sector comprised between those radii which appertain to the extreme 
points of the series. And when the circuit of the curve is wholly 
included in the summation, the radius PT having then passed round 
the circumference of the circle and again reached its original position, 
the total sum of all the segments, cut off by the sides of the polygon, 
will be practically represented by one-sixth of the entire circle so 
described. 

A higher degree of accuracy will obtain if the values thus ex- 
hibited be each of them increased by a small correction found by 

3 
dividing it by the number - 76n 2 — 1, or -n 2 — 1 nearly. 
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From the formula (13) and (14) it may be inferred that if the total 
length of the curve be divided into n equal portions, and a circle be 
described having the length of one of these portions as radius, then 
the total sum of the segments will be absolutely less in value than 
one-sixth of this circle, and greater than one-sixth of the circle 
whose radius is equal to the side of the polygon. 

Note. — The general development given in the foot note of page 317 
may be obtained by the process there indicated. The following is an inde- 
pendent investigation of it by means of the differential calculus. 

To determine the value of the curvilinear segment which connects two 
consecutive ordinates V and V, , it will be convenient in the first instance 
to take the origin midway between them; then, by Maclaurin's theorem, 

„ Tr ,{dV\ , /ri 2 V\x 2 fd s V\x 3 (d*V\ x* „ 

in which the differential coefficients also apppertain to the middle point. If 
Sx denote the intervening distance between consecutive ordinates, the values of 
V and V, will be obtained by respectively putting — \lx and -|--g&£ for x. 
Hence, taking the difference and half sum of the values thus found, we get 

„~ ,r ,(dW\Sx* /J 4 V\ Sx* (d 6 V\ Sx* „ 

And as the formula (19) will apply generally to any other function, it will 
give, for the curvilinear area contained between V and Vj , 

wrojl „. /dw\ gx3 /d*v\ a* 5 fd 6 v\ Sx? 

S{/ Ydx ) =y iS x + ^)— i+ {—) ¥: ^+{-^) ¥ --^ . (21). 

Also the trapezoidal area between V and V! is 

,™» ^ * , fd^V* 3 , fd*Y\ Bx 6 /d 6 V\ W „ , nn . 

(V)S*=V^ + (— j_ +(_ ) §XIi - 2 + (— ) 5-^ (22). 

The curvilinear segment is found by subtracting (22) from (21), and is 

Segment=- [— )— - {^^j, _^_ j___ fa. . (23). 

But, according to the formula (19), 

.fdV\ ,dW\. ,'d*V\ Sx* idW\ lx? 

\-dj)= {ur + {ix«)2jA +&c - 

•(£)- o- +*■ 

&c. &c. &c. 
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Therefore, by substitution in (23), 

Se gm ent = _-S(-j + 3 - ixr ^(^j- ——g^ — ] & c. .(24). 

and as the interval dx is constant, we hence obtain, by summation from 

2(segments) ....„= _)(--) _^--)J 

3.4.5.12 1 V «k" /o \ ^ /.) 



V to V„ 



3. 



+ 3 



...7.12 I Ui 8 J V «&•/.' 
&c. &c (25). 

If the interpolated interval — be substituted for fix, that is, if the 

respective terms of (25) be divided by m 2 , m*, m\ &c, it will then 
express the value of the sum of the interpolated segments between 
the same limits. Now, as the total curvilinear area is made np either 
of the polygon 2(V)<5.r and the segments (25), or of the interpolated 

polygon 2(V') — and the corresponding interpolated segments, it will 

follow that 

S(V')— =S(V)S*+S(segments)— 2 (interpolated segments) 
m 

■ 2( £— oo + £H£H£).} 

'20m* \{dx 3 ) \dx*)„\ 



m 

TO' — ! 



720m 4 
m 6 -l 



*<m-m 



1 30240m 6 

&c. &c (26). 

In the corresponding formula on page 317 the value of Sx is supposed 
to be unity. 

(To be continued.) 



